An Electron Spin Resonance Selection Rule for Spin-Gapped Systems 
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The direct electron spin resonance (ESR) absorption between a singlet ground state and the 
triplet excited states of spin gap systems is investigated. Such an absorption, which is forbidden 
by the conservation of the total spin quantum number in isotropic Hamiltonians, is allowed by the 
Dzyaloshinskii-Moriya interaction. We show a selection rule in the presence of this interaction, 
using the exact numerical diagonalization of the finite cluster of the quasi-one-dimensional bond- 
alternating spin system. The selection rule is also modified into a suitable form in order to interpret 
recent experimental results on CuGeOs and NaV2 05. 
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The spin gap generated by strong quantum fluctua- 
tions in low-dimensional systems, such as the bond alter- 
nating chain, the spin 1-chain or spin ladders, continues 
to attract much interest, both experimental and theo- 
retical. Electron spin resonance (ESR) experiments at 
very low temperatures provide very high-energy resolu- 
tion measurements of the spin gap. Tuning an external 
magnetic field, the field-dependent energy of the excited 
triplet state is adjusted to the frequency of the electro- 
magnetic propagating wave. A direct transition from the 
singlet ground state to the gapped triplet state may be 
observed. Such a direct transition would be, however, 
forbidden by the conservation law of the total spin quan- 
tum number for a perfectly isotropic but gapped spin 
system. 

In recent ESR measurements for the spin-Peierls sys- 
tem CuGeOs O and the bond-alternating chain system 
NaV205 0, direct transitions from the ground state to 
the first excited triplet state have been detected. Since 
the observed intensity of the transition for CuGeOa was 
almost independent of the external magnetic field, Nojiri 
et al. ||I| concluded that the explanation may be provided 
by the Dzyaloshinskii-Moriya (DM) interaction, which 
had been suggested on the basis of electron paramagnetic 
resonance measurements |3||. In contrast, another mech- 
anism based on the effective staggered field due to alter- 
nating g-tensors Q , which had successfully explained the 
ESR absorption of the Haldane gap for the S = 1 antifer- 
romagnetic chain Ni(C2lIgN2)2N02(C104), abbreviated 
NENP, should lead to a strong field dependence of the 
intensity. In addition, Kokado and Suzuki's phenomcno- 
logical Hamiltonian including the DM interaction repro- 
duced some of the results from the angle-dependent ESR 
measurement of CuGeOs 0. These works motivated the 
derivation of a selection rule for direct ESR absorption in 
the presence of the DM interaction. In the present letter, 
we present such a general rule which emphasizes the role 
of the relative orientations of the external field, the po- 
larized propagating magnetic field and the vector of the 
DM interaction. We also present a more suitable form for 



the experimental results, for which the electromagnetic 
wave is not polarized. 

When the magnetic field of the wave h(t) is polarized 
along an axis a {a = x,y,z), the intensity of ESR for 
quantum spin systems at zero temperature, assumed due 
to the magnetic dipole transitions, is determined by the 
matrix element 
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where 5" is the a-component of the local spin operator 
at the j-th site, \g) and |e) denote the ground state and 
the excited state with the energy difference E^ — Eg corre- 
sponding to the frequency uj of the electromagnetic wave. 
For an isotropic spin system, the dynamics conserves the 
total spin quantum number ^totai- Therefore, every com- 
ponent of the intensity (H) vanishes. In the presence of 
the DM interaction, however, S'totai is no longer a good 
quantum number and some intensities may become fi- 
nite. To investigate the effect of the DM interaction on 
the intensity, we first consider the problem of two coupled 
S = 1/2 spins with the Hamiltonian 
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JSi-S2 + D-(Si xS2)-H-(Si + S2 
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where D is the vector of the DM interaction and H is 
the external magnetic field. Solving the problem of two 
spins, we calculate the intensities of the ESR transitions 
from the ground state to the excited states. For H = 
and D II z, they are given by: 
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For nonzero magnetic field, this leads to the following 
selection rule. Nonzero intensity appears only when the 
h(i) has some component satisfying one of the following 
two conditions: 



(i) h(t) _L H II D: the intensity is independent of H . 
(ii) h(t) II H ± D: the intensities depend on H = |H| 
and the intensity of the transition to the nondegenerate 
state (iS'^otai = at iJ = 0), which vanishes aX H — Q 
increases as iJ^ for small H . For the degenerate states 
(S'j^Qjj^i = ±1 at i/ = 0), the intensities decrease. 

On the other hand, the intensity would still vanish if 
h(t) II H II D because of the rotational symmetry around 
these vectors. Since these rules depend only on the rel- 
ative configuration of the vectors h(i), H and D, and 
they are independent of the lattice structure, they are 
expected to be valid even in general bulk systems with a 
spin gap. (The if-independent intensity of rule (i) is gen- 
erally justified by the conservation of ^totai-) We will jus- 
tify them with the following numerical analysis on finite 
clusters of the quasi-one-dimensional bond-alternating 
spin system. In particular cases, special symmetries may 
lead to vanishing intensities at certain frequencies, in 
which case the gap observed in ESR may not coincide 
with that of neutron scattering. 

In order to examine the above selection rules in many- 
spin systems, we consider a two-dimensional lattice of 
weakly coupled bond- alternating chains. CuGeOa is one 
example of such a lattice. The cluster is shown in Fig. 
1, where the axes parallel and perpendicular to the chain 
correspond to the c- and ^axes of CuGeOs, respectively. 
The Hamiltonian with an external magnetic field is given 

by 
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FIG. 1. Cluster of 16 spins with DM interaction, used in 
the exact numerical diagonalization. The lattice is taken to 
be similar to that of CuGeOa . 



where the index I specifies a chain. The parameters are 
fixed as J=l, S = 0.1 and Jl=0.15, realistic for CuGeOa 
1^ . We consider the interchain and intrachain DM inter- 
actions (the direction of the vectors D is denoted z): 



"^DM.i — ^D^/j. (S/j X S; + ij) 
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Using the numerical diagonalization based on the 
Lanczos algorithm, we obtain the eigenvalues and eigen- 
vectors of the singlet ground state and the first excited 
triplet of the finite cluster described by the Hamiltonian 
(H) including the DM term. Under an external magnetic 
field, the triplet splits into three modes with energies 
A+, Aq and A_. To estimate the ESR intensity for 
each singlet-triplet transition, we calculate the matrix 
elements (||) for each excitation mode, denoted by /", 
Iq and /". The numerical calculation is performed for 
the 4 X 4-spins cluster shown in Fig. 1 with periodic 
boundary conditions in both directions. 



At first, the interchain DM interaction (g) is in- 
vestigated. Since the particular pattern of D-vectors, 
D\^i = (— 1)'D5_, leads to zero intensity for every polar- 
ization of the wave, we consider the sign-alternating one 
-1)'+^!?^. In the presence of the DM interac- 



tion with _D^ = 0.1, the ESR intensity is calculated under 
various external fields. To examine the selection rule (i), 
all the nonzero intensities are plotted versus the external 
field H^- which is parallel to D_l in Fig. 2(b). Figure 2(a) 
shows the i?^ dependence of the three modes which still 



have the good quantum number S^^^^y. Note that the 
triplet is split even in zero field because of the DM inter- 
action. Figure 2(b) shows that a finite intensity appears 
only when the polarized field has a magnetic component 
perpendicular to the DM vector and the external field. 
Since the magnetization along the external field is con- 
served in this case, the intensities are independent of iJ^, 
as shown in Fig. 2(b). These results agree with rule (i). 
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FIG. 2. (a) Energy gaps and (b) ESR intensities plotted 
versus the external magnetic field parallel to the interchain 
DM vector. 



The D^ dependence of the intensity is shown in Fig. 
3(a). The ln{Il)-\n{Dl) plot in Fig. 3(b) is almost hnear 
with a slope of 2. It indicates the relation 



P - \D 
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(8) 



which is satisfied in the two-spin system discussed above. 
(We omitted the symbol _L for D\^ in eq. (||), because it 
is also held for the intrachain DM interaction.) 
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FIG. 3. (a) ESR intensity I^ = II = II = P_ versus the 
amplitude of the DM vector 73^ (b) Plot of ln(/^)-In(L»^). 
It suggests 7+ ~ [D'^f. 



The results for the external field {H^) perpendicular to 
the DM vector are shown in Figs. 4(a) and 4(b). Figure 
4(b) shows that the intensities at the energies A4. or A_ 



are finite only if the polarized field h(t) has one compo- 
nent along the external field. The intensity /^ increases 
as H^ increases. These behaviors support the validity of 
rule (ii) even in many-spin systems. 
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FIG. 4. (a) Energy gaps and (b) ESR intensities plot- 
ted versus the external magnetic field perpendicular to the 
interchain DM vector. 



Next, the DM coupling along the chain (0) is investi- 
gated. The intrachain DM interaction gives rise to some 
nonzero intensities in the case of a uniform sign, that is, 
PI 



IKj 



= PI. The results for such an intrachain DM in- 



teraction with PI =0.1 are presented in Figs. 5 and 6, 
for the cases of H || D and H _L D, respectively. These 
results reveal the same behaviors as the interchain DM 
interaction for all the wave polarizations, although the 
absolute value of the intensity is slightly larger. Thus, 
selection rules (i) and (ii) are revealed to be valid also 
for the DM coupling along the chain. These results sug- 
gests that the rules do not depend on where the DM 
interaction acts. 
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FIG. 5. (a) Energy gaps and (b) ESR intensities plotted 
versus the external magnetic field parallel to the intrachain 
DM vector. 
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FIG. 6. (a) Energy gaps and (b) ESR intensities plot- 
ted versus the external magnetic field perpendicular to the 
intrachain DM vector. 



We derived general selection rules (i) and (ii) for a 
polarized electromagnetic wave. In the experiment on 
CuGcOa, however, the wave was not polarized, and its 
wave vector k was oriented along different crystal axes. 
In order to compare with these experimental results, we 
would average the calculated ESR intensities over the 
two components of the magnetic field perpendicular to 
the wave vector k. In the following, we discuss the mod- 
ified rules for nonpolarized electromagnetic waves, in the 
Faraday or Voigt configurations. 

(A) Faraday configuration: k || H 

Since the magnetic field h(i) perpendicular to k is always 



perpendicular to the external field, ESR absorption is al- 
lowed only if the rule (i) is satisfied. The intensity should 
then be independent of the external field. It is expected 
to be largest for H || D and zero for H ± D. Based on 
the eq. (H), the measured intensities for the cases where 
the external field is successively along a-, h- and c-axes, 
would enable us to determine the direction of the vector 
D in the following way. The intensity at the energy A-|_ 
for H ||a in the Faraday configuration is noted l%\-, so 
that the components of D should be determined by 
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(B) Voigt configuration: k ± H 

Since h(i), which is perpendicular to k, has components 
parallel and perpendicular to H, both rules (i) and (ii) 
can be realized. We discuss the two possibilities sepa- 
rately. 

(a)H II D (rule (i)): every intensity is independent of 
the external field and the intensities at the energies A+ 
and A_ are half that of the Faraday configuration with 
H II D. 

(/9)H ± D (rule (ii)): the ESR intensities at the energies 
A+ and A_ depend on the external field. 

In a vanishing external magnetic field, symmetry argu- 
ments based on the rotation of the spin operators H] can 
be given. In particular, we explain the vanishing inten- 
sity found for a transverse alternating DM interaction as 
well as the D^ dependence of the intensity for the others. 

First, consider the alternating transverse DM inter- 
action D^i = (—lyD^. Following Oshikawa and Af- 



fleck 



we introduce rotated spin operators S- 
Sf, with ia.n6 = D 



S/2c± 



QlZ 



l± _ 
J 

This 



exactly maps the Hamiltonian onto an XX Z model with 
alternating coupling for S' operators for which we can 
neglect the small second-order anisotropy. Transforming 
the Zeeman coupling gives rise to a staggered magnetic 
field h**°9ff(i) for rotated operators §. The intensity of 
the transition is given for the rotated variables as 
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The first term, involving the total spin operator, still 
gives a vanishing contribution because the rotated Hamil- 
tonian is isotropic. The second term comes from the 
field that is staggered in the transverse direction. It 
therefore has q± = tt symmetry. On the other hand, 
the lowest-lying excitation and the ground state have 
q± — symmetry for an antiferromagnetic transverse 
interaction for the lattice of Fig. H. |q| This is why the 
transition from the ground state to the first gapped ex- 
cited state is still forbidden with this DM anisotropy. 
A transition to the triplet excited state q± = tt with 
higher energy is, however, allowed with an intensity fac- 
tor sin2(6l/2) = {Dl/2J±y in the small D limit, g 



Second, consider the case of a constant DM interac- 
tion Df.i ■ = Df.. A rotation to the new set of operators 
gi± ^ e^i^JS-i ^ s'l^'j = Sl^ (with tan 6* = Dy J ) also 
maps the Hamiltonian onto a bond-alternating XXZ 
model. In the rotated frame, the magnetic field h(i), 
however, transforms into a field which rotates in space 
with a periodicity Itx jd. The intensity is then: 
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Therefore, the intensity of the ESR transition at g = 
induced by uniform DM interaction equals the intensity 
at small q = 6 calculated for a model without DM inter- 
action. For small g, the intensity should behave as / ~ g^ 
(the matrix element should be analytic for a spin-gapped 
system), so that /^'^ ^ (P = (DljjY, in agreement with 
the relation found numerically. 

Finally, we briefly discuss the experimental results of 
CuGeOa on the basis of the above selection rules. We as- 
sume that the forbidden absorption observed is induced 
by a DM interaction and see what are the consequences. 
Considering only the absorption at the energy A+, the 
measurement in the Faraday configuration revealed that 
Jfp-j is the largest one. Relation (0) thus indicates that 
D is along the c-axis. The interchain DM interaction 
with Dj^ along the c-axis is in fact the most realistic sit- 
uation with regards to the crystallographic structure of 
CuGeOs. The sign configuration D^^^^ — {—\)^D\ is, 
however, expected [^0, but it leads to zero intensity 
for every polarization, as proven by the symmetry argu- 
ment, and cannot explain the transition observed. The 
other patterns for interchain D\i = {—Xf^'-'D\ or in- 
trachain Dyj . = D\^ give intensities (the observed inten- 
sities would then indicate an extra small component of 
D along the a-axis), but we cannot conclude from these 
measurements alone whether they exist or whether we 
need other interactions such as exchange anisotropics or 
alternation of the (7-tensor to explain the ESR absorp- 
tions. In particular, the measured intensities in the Voigt 
configuration cannot be compared with the above rules, 
because the external field dependence has not been mea- 
sured yet. The observation of increasing intensity with 
the external field would be good evidence of the contri- 
bution of the DM interaction to the absorption. Here, 
we can make the following observation for the case of 
NaV205: An absorption at the energy gap with no strong 
field dependence has been observed in the Faraday con- 
figuration H for an external field along the a-axis. If we 



assume that it is mainly due to DM interaction, the above 
rule (A) would tell us that the DM vector is along the 
a-axis, in agreement with the crystallographic structure 
of V-O-V bonds along the chains. 

In summary, we presented two selection rules (i) and 
(ii) for the ESR polarized absorption for a spin-gap sys- 
tem in the presence of the DM interaction. For certain 
frequencies there may be additional selection rules, de- 
pending on the spatial variations of the DM interactions. 
We also modified the rules for a nonpolarized electro- 
magnetic wave for comparison with recent experiments. 
It was found that the intensity depends on the external 
field in some of the Voigt configurations. 
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